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ABSTRACT 

We present an estimate of the global budget of baryons in all states, with 
conservative estimates of the uncertainties, based on all relevant information we have 
been able to marshal. Most of the baryons today are still in the form of ionized gas, 
which contributes a mean density uncertain by a factor of about four. Stars and their 
remnants are a relatively minor component, comprising for our best-guess plasma 
density only about 17% of the baryons, while populations contributing most of the blue 
starlight comprise less than 5%. The formation of galaxies and of stars within them 
appears to be a globally inefficient process. The sum over our budget, expressed as a 
fraction of the critical Einstein-de Sitter density, is in the range 0.007 <J £Ib ^ 0.041, 
with a best guess CIb ~ 0.021 (at Hubble constant 70 km s _1 Mpc" 1 ). The central 
value agrees with the prediction from the theory of light element production and 
with measures of the density of intergalactic plasma at redshift z ~ 3. This apparent 
concordance suggests we may be close to a complete survey of the major states of the 
baryons. 



Subject headings: cosmology: observations — galaxies: fundamental parameters 



1. Introduction 

The evolution of structure in the expanding universe has redistributed the baryons from 
a nearly smooth plasma at the time of light element nucleosynthesis to a variety of states — - 
condensed, atomic, molecular, and plasma — in clouds of gas and dust, planets, stars, and stellar 
remnants, that are arranged in the galaxies, in groups and clusters of galaxies, and in between. 
The amounts of baryons in each state and form at low redshift can be compared to what is 
observed at higher redshift and to the total cosmic abundance predicted by the theory of element 
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production in the early universe. Budget estimates must be informed by ideas on how structure 
evolves as well as by the observations, and the budget in turn is a test of these ideas. Knowledge 
of the baryon budget is an essential boundary condition for the analysis of how structure formed 
and of the nature of cosmic dark matter. 

Advances in observations allow reasonably sound estimates of the amount of baryons present 
in a considerable variety of forms. In a previous discussion (Fukugita, Hogan, & Peebles 1996), 
we presented a picture for cosmic evolution suggested by the results of budget estimates. In this 
paper we give details of the budget calculations and update them using our choices for the current 
best knowledge relevant to the calculation. Persic & Salucci (1992), Gnedin & Ostriker (1992), 
and Bristow & Phillipps (1994) estimate baryon abundances with different emphases; the results 
are compared in §5.1. 

The main focus of this paper is the baryon budget at low redshift; our accounting is presented 
in § ^. The low redshift budget can be compared to the situation at z ~ 3, where quasar absorption 
lines allow a comprehensive accounting of the diffuse components. We comment on this in § ||. The 
budgets are compared to the constraint from light element production in § [|, and the implications 
for galaxy and structure formation are summarized in § [5[ 

We write Hubble's constant as 

H = lOO/i km s _1 Mpc -1 . (1) 

Where not explicitly written we use solar units for mass and luminosity and megaparsecs for the 
length unit. 

2. The Baryon Budget at z ps 

2.1. Stars and Remnants in Galaxies 

Stars in high surface density galaxies are the most prominent place for baryons. For our 
purpose it is reasonable to imagine the galaxies contain two distinct stellar populations, an old 
spheroidal component and a disk component consisting of generally younger stars, with a mix 
depending on galaxy type. Elliptical galaxies lack a significant disk component, and irregular (Im) 
galaxies are in the opposite extreme, having small or insignificant bulges. The formation of these 
two components seems to follow different histories, so it is appropriate to count baryons in stars 
divided into these two categories. The mass density for each component is obtained as 

/)(sph, disk) = L B ■ /s(sph, disk) • (M/L B } sp h,disk, (2) 

where £ is the mean luminosity density, fg is the fraction of the luminosity density produced 
by the spheroid or disk component, and {M/Lb) is the mass-to-light ratio for each component 
including the stars and star remnants. The suffix B refers to luminosities measured in the B band, 
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our choice for a standard wavelength band. Although the stellar population of irregular galaxies 
is similar to that of the disk component, we treat the former separately to emphasize its distinct 
role in the luminosity density at low redshift and a slightly smaller mass-to-light ratio because of 
its younger mean age. 

2.1.1. Luminosity Fractions 

We consider first the luminosity fractions / in equation (^). In Table |l| we show the bulge 
fraction of galaxies in luminosity, k = -^bulge/(-^bulge + Aiisk)i f° r each morphological type, from 
the bulge-disk decomposition of Kent (1985) based on his accurate Gunn r band CCD photometry. 
We take the median of the distribution presented by Kent. The bulge fraction in r band is 
converted into that in the B band using the color transformation law B — r = 1.19 for elliptical 
galaxies and bulges (Fukugita, Shimasaku & Ichikawa 1995). The values of n thus obtained in 
the B band are consistent with the result of Whitmore and Kirshner (1981), although the latter 
sample is not as large as that of Kent and the decomposition, based on photography, may be 
less accurate. The bulge fraction is well correlated with the morphology (Morgan 1962; Dressier 
1980; Meisels & Ostriker 1984; Kent 1985), but the morphologies at given bulge fraction overlap 
neighboring classes. 

Row (3) in Table | is our adopted distribution of morphologies, expressed as the fraction 
\it of the mean luminosity density contributed by galaxies of morphological type t. Postman 
&; Geller (1985) give E:S0:S+Irr=0. 12:0.23:0.65. This ratio is consistent with the Tinsley mix 
(Tinsley 1980), which subdivides the spiral fraction into further subclasses, while E and SO are not 
separated. The luminosity fractions in equation (|2|) are obtained by combining the two results: 

/ B (sph) = £> tKt = 0.385, / B (disk) = X> t (l-K t ) =0.549, / B (irr) = /x irr = 0.061. (3) 

The sum in the second expression is from SOs through spirals. Because the mass in SO discs is so 
small we make little error in treating the disks of spirals and SOs as the same population. The 
separation into each morphological type is ambiguous, and the classification depends much on the 
author. To model such error, we recalculate the fractions / after reclassifying 1/3 of the galaxies 
in each type systematically to the neighboring later type, or else to the neighboring earlier type. 
This produces the range (0.324~0.459):(0.599-0.499): (0.076-0.041). When we use a morphological 
fraction or a bulge-to-disk ratio taken from other authors, we find the luminosity fractions / fall 
within this range in most cases. In particular these estimates, when translated into the V pass 
band, agree with the survey of Schechter and Dressier (1987). This spread of values of the fit is 
taken into account in our assignment of errors in estimates of the mean mass density. 
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2.1.2. The Mean Luminosity Density 

Recent measurements of the mean luminosity density contributed by high surface brightness 
galaxies at moderately low redshift are summarized in Table ^. There is a substantial spread of 
estimates of the parameters <f>* and M* in the Schechter function, but the total luminosity density, 
C = (j>*L*T(a + 2), is better defined: most values are in the range Cb = (1.8 — 2.2) x 10 8 /iL Q Mpc~ 3 . 
The exception is the significantly lower number from the APM-Stromlo survey (Loveday et al. 
1992). We adopt 

C B = (2.0 ± 0.2) x 10 8 /iL o Mpc" 3 , (M/L B ) crit = (1390 ± U0)h. (4) 

The second expression is the corresponding value of the mass-to-light ratio in the Einstein-de Sitter 
model. 

2.1.3. Mass-to-light Ratio of the Spheroid Component 

The mass-to-light ratio of stars and their remnants in spheroid populations can be inferred 
from star population synthesis calculations, assuming the stars are close to coeval (so that the 
integrated B luminosity is not much affected by recent episodes of star formation) and the 
spheroid stars in our neighborhood are representative. Within these assumptions the most serious 
uncertainty in the calculation is the lower cut off of the initial mass function (the IMF). In the 
past the cutoff was somewhat arbitrary, between O.IMq and O.2M0, and the choice seriously 
affected the estimate of M/L. Thus for the Salpeter IMF M differs by 24% between these two 
choices, while L is virtually unaffected. This uncertainty for the local IMF has been removed by 
the recent advance in observation of M subdwarf star counts with the Hubble Space Telescope: 
it has been shown that the IMF of the local Galactic stars has a turnover at about O.4-M0 
(Gould, Bahcall & Flynn 1996; hereafter GBF), and the integral over the GBF mass function 
converges without an arbitrary cut off. For M > M & the GBF mass function agrees well with 
the Salpeter function. Using the GBF mass function for M < Mq, instead of the Salpeter mass 
function with x = 1.35 (where x is the index of the IMF dn/dM oc M~^ +x '), the mass integral 
is 0.70 times that obtained with the Salpeter mass function cut off at 0.15M@. The light in the 
blue-visual bands contributed by stars less massive than one solar mass is small: from Table 4 
of Tinsley & Gunn (1976), we estimate that the V light integral decreases only 2% with the use 
of the GBF mass function. Therefore, the mass-to-light ratio with the GBF mass function is 
(M/Lv)gbf = 0.72(M/Ly)g a i peter , where the latter ratio is computed for the cutoff at O.15M . 
Chariot, Worthey, & Bressan (1996) have given a compilation of mass-to-light ratios from various 
stellar population synthesis calculations. They show that MjLy is reasonably consistent among 
authors when the lower mass cutoff is fixed. From the Chariot et al. (1996) compilation of 
population synthesis calculations, the GBF IMF, and the conversion formula discussed above, we 
obtain (M/Lv)gbf = (4.0 ± 0.3) + 0.38(tc — lOGyr), where to is the age of the spheroid and 
the error represents the model-dependent uncertainty. Using the Chariot et al. (1996) synthetic 
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calculation of the B — V color, we obtain M/L B = (5.4 ± 0.3) + 0.7(t G - lOGyr). For t G = 12 ± 2 
Gyr we have M/Ly = 4.0 - 5.9, or 

M/Lb = 5.4 — 8.3. (5) 

The M/L ratio in equation (]5|) may be compared to the values inferred from the kinematics 
of the nuclear regions of elliptical galaxies, van der Marel (1991) finds M/Lr = (6.64 ± 0.28)/i 
(Johnson R), which translates to M/Ly = (7.72 ± 0.33)/i or 

M/Lb = (9.93 ± 0.42)/i, (6) 

using (V — Rj) = 0.68, corresponding to the average color (B — V) = 0.92 for the 37 elliptical 
galaxies he used ((B — V)q = 0.65 and (V — Rj)q = 0.52 are also used). There is excellent 
consistency between equations (|5|) and (||) if the Hubble constant is in the range h = 0.6 — 0.8 and 
spheroid populations dominate the the mass of these nuclear regions. 

Dynamical measures show M/L depends on the luminosity in elliptical galaxies (Kormendy 
1986): van der Marel (1991) finds M/L = (L/L*) a35±a05 . This manifestation of the color- 
magnitude relation of early-type galaxies probably reflects the effect of metallicity and perhaps 
also of dark matter halos. In the former case we estimate this dependence may reduce the 
effective mean value of M/L weighted by the contribution to the luminosity density by 17%, and 
accordingly reduce the lower end of the allowed range by this amount. 

We conclude that a good estimate of the mass-to-light ratio of stars and their remnants in 
spheroids is 

M/Lb (spheroids) = 6.5 + ^. (7) 



2.1.4- Mass-to-Light Ratio of the Disk Component 

One way to estimate the mass-to-light ratio of the disk stars is to use estimates of the mass 
and luminosity column densities from surveys of stars and remnants in the solar neighborhood. 
Commonly used values are T,m — 45M pc -2 (stars and stellar remnants) and ~ 15L Q pc~ 2 
(Binney & Tremaine 1987; Bahcall & Soneira 1980). GBF find that the mass column density must 
be reduced to S^,/ ~ 27.3M Q pc~ 2 because of the turnover of the luminosity function towards the 
faint end, as noted above. With this new value we have M/Ly ~ 1.82 or M/Lb — 1-50 for the 
disk with B — V ~ 0.44 (de Vaucouleurs & Pence 1978). This is not very different from the early 
estimate by Mihalas & Binney (1981), M/L B ^ 1.2. 

The solar neighborhood value may be compared with M/Lb from population synthesis 
calculations. The Tinsley (1981) 10-Gyr model with a constant star formation rate gives 
M/Lb = 1.2 after a 20% upward correction in the mass integral of the IMF to match GBF. 
Shimasaku & Fukugita's (1997) gas infall model (with effective disk age ~ 10 Gyr) gives 
M/Lb = 1.9 after a 40% reduction in M to correct the IMF. From a study of rotation curves for 
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spiral galaxies, Persic and Salucci (1992) find M/L^ = 1.24/t. We conclude that there is reasonably 
small scatter among these different approaches to the disk mass-todight ratio, and we adopt 

(M/L B ) disk = 1.5 ±0.4. (8) 

There is good evidence that stars in gas-rich irregular (Im) galaxies are young. (For reviews 
see Fukugita, Hogan, & Peebles 1996; Ellis 1997). If the mean age of irregular galaxy stars is 
~ 5 Gyr, the B luminosity is 0.5 mag brighter and the mass in star remnants is about 10% smaller 
than for a 10 Gyr age disk. These corrections reduce M/L in equation (@) to M/Lb = 1.1 for 
irregular galaxies. A population synthesis model for 8 Gyr age (Shimasaku & Fukugita 1997) also 
gives M/Lb = 1-1 after the correction to the IMF. We take 

{M/L B )\n = 1-1 ±0.25. (9) 



2.1.5. Baryons in Stars and Remnants 

With the above numbers the mean mass densities in our three classes of stars and their 
associated remnants are 

/ ±0.00121 \ , . . 

^spheroid stars = f 0.00180 _ Q Q0Q ) h' 1 , (10) 

ttdisk stars = (0.00060 *™°°2°) ^\ (11) 
^stars in Irr = (0.000048 +0 - 00033 ) ( 12 ) 

V -0.00026/ ' v ; 

in units of the critical density p cr it = 2.775 x lQ ll h 2 M & Mpc -3 . Here and below the upper 
(lower) errors are obtained by taking all errors in each step to maximize (minimize) the result. 
In this sense the error ranges are conservative. The numbers in equations ([h]) to ( |i~2|) depend 
through the M/L calculation on age estimates, and hence depend on the Hubble constant h and 
other cosmological parameters in a complicated way. If dynamics were used to estimate M/L the 
estimates of Vt would not depend on h, but in that case h ~ 0.7 would be needed for dynamics to 
agree with the synthesis calculations. Either way, consistency holds for h ~ 0.7 in a low density 
universe. 



2.2. Atomic and Molecular Gas 

From HI 21-cm surveys Rao and Briggs (1993) find Hhi = (1-88 ± 0.45) x 10~ 4 /i _1 (see also 
Fall & Pei 1993, Zwaan et al. 1998). A 32% addition (for 24% helium mass fraction) to the 
Rao-Briggs value for helium yields 

^atomic = (0.00025 ± 0.00006)/i _1 . (13) 
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This result is quite secure because absorption line studies show directly that most of the HI is in 
systems with column densities high enough to detect in 21-cm emission. These columns are also 
generally high enough to resist photoionization so we are counting here mostly regions dominated 
by neutral atomic gas. 

The ratios H2/HI of mass in molecular and atomic hydrogen listed in Table [l] for each 
morphological type are median values taken from the compilation of Young & Scoville (1991), 
who base the molecular hydrogen masses on CO observations. The representative masses in 
atomic hydrogen in Table |] are from Roberts & Haynes (1994). Combining these numbers with 
the distribution m in morphology, weighted by HI mass, we obtain the global ratio in galaxies 
PH 2 /phi = 0.81, and from equation (|l3|) 

n R2 = (0.00020 ± o.ooooe)/^ 1 . (14) 



2.3. Baryons in Clusters of Galaxies 

The cluster mass function is approximated as (Bahcall & Cen 1993) 

n cl (> M) = 4 x lO" 5 /^— ) exp ( - — ) Mpc" 3 , (15) 

where M* = (1.8 ± 0.3) x 1O 14 /i _1 M , and M is the total gravitational mass within a sphere of 
a radius of 1.5/i -1 Mpc (the Abell radius) centered on the cluster. At the Abell radius the mass 
distribution is close to dynamical equilibrium. The envelope of matter around the cluster at 
greater distances merges into the "large-scale structure" and we take this matter to be included 
in the field (as discussed in §2.4). We define clusters as objects with mass M > 10 14 hM & . The 
integral J dM Mdn c \j 'dM gives the mean mass density in clusters, 

Pd = ( 7 - 7 ^_ 2 2 l) x W 9 h 2 M Q Mpc" 3 . (16) 

The contribution of this gravitational mass to the density parameter is 

tt*- 0.028 +a °° 9 . (17) 
-0.008 v 1 



Intracluster plasma masses are well determined from X-ray observations (Fabricant et al. 
1986; Hughes 1989; White et al. 1993). The ratio of X-ray emitting gas mass to gravitational 
mass within the Abell radius from the survey by White & Fabian (1995) is 

(M Hn /M grav ) cl = (0.056 ± 0.014)/T 3 / 2 . (18) 

In the White & Fabian (1995) sample this ratio shows no correlation with cluster mass. The value 
(18) is independently verified by Myers et al. (1997), (M H n/M grav ) c i = (0.061 ± O.Oll)/^ 1 from 
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the measurement of the Sunyaev-Zeldovich effect in three clusters. The product of equations (|1 
and (jlj) is 

(^n) cl =fo.00155 +a001 °°^--. (19) 
v miyci ^ -0.00072/ v ; 

This is the contribution to the baryon budget by intracluster plasma. 

Let us estimate the baryons in stars in galaxies in clusters. Although these stars are included 



in the estimate for stars in § 2.1, it is useful to compare their mass with the plasma mass, which 
we can derive from the cluster mass-to-light ratio. In the discussion of the Coma cluster by White 
et al. (1993), straightforward applications of galaxy velocity dispersions or the X-ray pressure 
gradient give (M gI . av /Ls)ci ~ 370/t, or values as large as 500/i if the analysis is constrained by 
models from numerical simulations of cluster formation. The CNOC value (Carlberg et al. 1996) 
transformed to the B band by (B — r)@ = 0.65 and {B — r) = 1.03 ± 0.1 for the average color 
of SO galaxies, after a passive evolutionary correction for early type galaxies, is somwhat larger, 
560/i. However, the CNOC luminosity density is correspondingly smaller than our adopted value 
(equation (4)), and it is the product of this with M/L that matters in deriving most global 
quantities; moreover we suspect that the product is more reliably estimated since in the end we 
are estimating the luminosity density in cluster galaxy stars, which does not depend on a mass 
estimate. We therefore adopt a central value, 

(M grav /L B ) cl = (450 ± 100)/i. (20) 

For cluster galaxies (Dressier 1980) we adopt the composition within the Abell radius 

/i(E) = 0.21, /x(SO) = 0.44, fi(S) = 0.35. (21) 

(This is somewhat richer in early-type galaxies than the estimate given by Schechter & Dressier 
(1987), but their sample extends beyond the Abell radius.) The mass-to-light ratios in equations (|7|) 
and (H) and the bulge-to-disk ratios in Table |l| give (M/L) star = 4.5 ±1. The ratio to equation ( ^0|) 
gives /iM stars /Mg rav = 0.010 ^q*oo4" product with equation ([l7| ) is the contribution to the 
density parameter by stars in clusters, 

(CWs)cl= (0-0003+^^) ^ > < 22 ) 

about ten percent of the total in equations (|l0|) to (FT 



2.4. Baryons in Groups 

Plasma associated with galaxies outside the great clusters makes a significant and still quite 
uncertain contribution to the baryon budget. The detected X-ray emission from plasma in groups 
is softer than for clusters, and softest for groups dominated by late-type galaxies. It is a reasonable 
presumption that the absence of detections often represents a lower virial temperature, rather 
than the absence of plasma. 
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We follow two approaches in estimating the mass in plasma associated with galaxies outside 
clusters. In the first, we start with a characteristic value for the ratio of plasma mass to 
gravitational mass in groups with X-ray detections. We apply this ratio to all galaxies, on the 
assumption that all field galaxies are likely to have inherited a similar plasma mass, to estimate 
the density in the warm phase (§ pAip . To this we add an estimate from quasar absorption of 
the plasma mass in low surface density cool clouds that certainly are not counted in the X-ray 
measurements (§ [2.4.2 ), The second approach (§ |2.4.3| ) assumes the plasma-to-gravitational mass 
ratio is the same in clusters and the field. 



2. 4-1- Warm Plasma in Groups 

Some surveys are available from low energy X-ray observations with ROSAT (Mulchaey et al. 
1996). The X-ray emission from 18 groups with total mass ranging from 1.2 to 8.3 x 10 13 /iMq 
corresponds to plasma mass fractions ranging from 0.004/i~ 3//2 (H97) to 0.09/i~ 3 / 2 (NGC 4261). 
The average is 

(MHii/M grav ) group = (0.022 ± 0.005)/iT 3/2 , (23) 

significantly below the number for clusters (eq. JT^j). In fact the baryon fraction shows a trend 
increasing with the group mass, approaching to the cluster value at the high mass end. This could 
be because groups are intrinsically poorer in plasma, or because much of the plasma is cooler and 
so escapes detection as an X-ray source. The latter is in line with the shallower gravitational 
potential wells in groups. The cool plasma clouds detected by Lyman-a resonance absorption 
(§ 2.4.2) similarly are not detected as X-ray sources. Thus the plasma identifiable from its discrete 



X-ray emission might be considered a lower limit to the net plasma associated with groups of 
galaxies. 



To convert equation (|23| ) into a mean baryon density we need the mean gravitational mass 
associated with field galaxies (which almost always are in groups). The following measures may be 
compared. First, we can extrapolate the Bahcall-Cen (1993) mass function (eq. |lf|), which they 
determined for M > 10 13 /iMq, to systems with the mass characteristic of galaxies, M ~ 10 12 HMq. 
The result of integrating this mass function from M = W 12 hM@ to M = 10 14 /iMq is 



= 0.12 ±0.02. (24) 

The cutoff is the characteristic mass of an L* galaxy, the minimum for a group. 

Second, we have dynamical mass measures from analyses of systems of galaxies on scales 
smaller than about 10/i -1 Mpc and outside the rich clusters. The survey of results of these 
analyses by Bahcall, Lubin & Dorman (1995) indicates M grav /L ~ (200 ^q )/!. This with the 
mean luminosity density in equation (|j) gives 

O = 0.14 + °- 10 . (25) 
-0.05 v ' 
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Third, we can use the luminosity density, scaling from M/L calibrated in the great 
clusters by taking account of the difference in luminosities from the difference in morphological 
mixes (see also Carlberg et al. 1997). The small scatter in the color-magnitude relation for 
ellipticals and the spheroid components of spirals suggests these stars formed early, so it is 
reasonable to assume that the ratio of spheroid luminosity to gravitational mass is the same 
in clusters and the field. It would follow that the cosmic mean mass-to-spheroid-light ratio is 
(M/Lb)cos = (-^grav/-^s)ci/(sph, cos)//(sph, cl) = (270 ± 60) h, where /(sph,cos) is the cosmic 
mean spheroid population luminosity fraction (eq. the value in clusters, /(sph,cl) = 0.64, 



follows from equation (21), and (M grav / Lb) c \ is from equation (20). With equation (f|) we 



have Q = 0.19 ^j'?!. This however gives the cosmic mean including clusters; for comparison to 
equation (p5| ) we ought to convert to the value outside rich clusters. The ratio of equation ( |i~6| ) to 



equation (20) is the mean luminosity density provided by the clusters, C c \ = 0.09£. The spheroid 



luminosity fraction in the field is 

/fieid(sph) = (£/ cos - C cl f cl )/(C - C cl ) = 0.94/ COS . (26) 

The assumption that the mass-to-spheroid-light ratio is universal thus indicates the density 
parameter in gravitational mass outside the Abell radii of the great clusters is 

n-a«™ (27 ) 

Despite the substantial uncertainties in each of these arguments we are encouraged by 
the consistency of equations (|24]), fl25|), and ( p7| ) to conclude that the density parameter in 
gravitational mass that clusters with galaxies on scales ^ 10ft, -1 Mpc is likely to be in the range 

« = 0.15+°;J°, (28) 

and that the spheroid-to-dark matter and baryon-to-dark matter ratios indeed do not vary 
widely between clusters and the field. As the evidence has indicated for some time (Peebles 
1986) and has been widely noted in recent years, this low density parameter universe offers a 
natural interpretation of a variety of observations. Recent examples include the age/distance scale 
relation, the abundance of cluster baryons (White et al. 1993), the growth rate of correlation 
functions (Peacock 1997), the growth rate of the cluster mass function (Bahcall, Fan & Cen 1997), 
and preliminary indications from the distant supernova Hubble diagram (Garnavich et al. 1998, 
Perlmutter et al. 1998). 



The product of equations (|23| ) and (|28| ) is the estimate of total plasma identified by X-ray 
emission in groups, 

(^HIl)group = (O- 003 ^^) k ' 3/2 ' (29) 

We remark that this estimate decreases only by 20% when we incorporate explicitly 
the trend that the baryon fraction increases as the group mass, (MHii/M grav ) groU p ~ 
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0.056/i~ 3 / 2 (M group /0.6 x 10 14 M o ) for M group < 0.6 x 1O 14 M . As we have noted, an alternative 
interpretation is that the trend is only apparent, a result of less efficient detection of plasma 
around cooler groups. 



2.4-2. Plasma in Cool Low Surface Density Clouds 

In addition to the hot plasma detected in X-rays, there is a considerable amount of plasma 
in cool, thermally stable, photoionized clouds with temperatures less than about 20,000 K. These 
cool low column density clouds are detected by the Lyman-a resonance absorption line in quasar 
spectra. They may be found as far as 5 Mpc from a galaxy, and generally seem to avoid both the 
voids and strong concentrations of galaxies (Shull, Stocke, & Penton 1996). This is consistent with 
the view that such clouds were left behind in the gravitational assembly of groups of galaxies; 
those that would have belonged to cluster members were incorporated and shock heated by the 
assembly of the clusters, joining the hot phase. In our budget for plasma in groups we thus use 
the sum of the mean baryon densities in these clouds and in cluster baryons detected by X-ray 
emission (eq. p£|). 



The estimate of the baryon abundance in such clouds is very sensitive to model details such 
as ionization and geometry. Let us assume as an illustrative model that the clouds are uniform 
density, isothermal spheres in photoionization equilibrium, and the size is of the order of 100 kpc, 
determined empirically by common absorption between two nearby lines of sight. Then the mass 
of the clouds that yield neutral hydrogen column density Nm are estimated, following Shull, 
Stocke & Penton (1996), to be 

M Lya = 7.1 x 10 8 Af o J%[T/{2 x 10 4 )] a375 A^ 4 /2 (i?/100kpc) 5/2 [(3 + a)/4.5]- 1/2 (30) 

where J-23 is the ionizing UV flux at the Lyman limit in units of 10 _23 erg s _1 cm _2 sr _1 Hz~ 1 , 
with a power index a, T is the gas temperature, N14 is the neutral hydrogen column 
density in units of 10 14 cm -2 , and R is a typical cloud size. The neutral hydrogen fraction 
x = n m /nu ^ 6 x 10" 5 JZli 2 [T/{2 x 10 4 )]-°- 375 iV 1 1 4 /2 (i?/100kpc)- 1 / 2 [(3 + aj/A.h] 1 ' 2 . With the 
detection frequency of the absorbers at z ~ 0, dN c /dz = 0o(vri? 2 )c 2 / Hq ~ 86, the mass density of 
the clouds is 

O ^ nnnoo /a^ 2 X IQ 4 )] ' 375 ^ WlOOkpc) 1 / 2 

^Lya = - 0.00234 — — —j- . 31 

Pcrit h[(3 + a)/Lb\ L l z 

This assumes spherical clouds. If the clouds are oblate the density decreases by the aspect ratio, 
which we assume to be between 1/5 and 1. 

Direct simulations of cloud formation and absorption spectra permit a much more detailed 
though model-dependent picture of the clouds. It has been demonstrated that these simulations 
give good fits to the observations at z « 3 (see below) but they do not not yet yield useful 
approximations to the situation at z = 0. Thus the direct estimates of the density parameter in 
this component remain quite uncertain. 



-12- 



2.4-3. Scaling from Clusters to Groups 

As an alternative to adding the previous two directly observed components (warm X-ray 
emitting gas and cool Lyman-a absorbing gas), we can estimate the total plasma in groups from 
the assumption that it has the same ratio to dark matter as in the rich clusters. That is, the 
scaling used to arrive at the density parameter in gravitational mass in the field from spheroid 
luminosity (eq. p7| ) readily generalizes to a scaling estimate of the density parameter in plasma 
associated with groups. The product of equation (|i~8| ) for the ratio of plasma to gravitational mass 
in clusters and equation (|28|) for the density parameter in the field yields 

(n H ii) fle id= fo.oos+^V 372 - (32) 



-0.004 

In a calculation along similar lines Carlberg et al. (1997) find a similar result, 
(^Hii)fieid = (0.01 — 0.014)/i~ 3 / 2 . The central value in equation ( |32[ ) agrees with the 
upper bound from the sum of the more direct estimates in equations (p9[) and pi]) . 



Equation ( p2[ ) assumes a fair sample of plasma relative to gravitational mass is assembled and 
preserved in the great clusters. If disks were less common in clusters because star formation has 
been suppressed, equation (]3^ ) would include baryons in disks in the field, but the correction is 
small because the baryon mass in disk stars is small. It is not thought that cluster cooling flows or 
winds have seriously depleted the cluster plasma mass. In a more direct argument Renzini (1997) 
finds that the iron and gas content of clusters and groups indicate the clusters indeed have close to 
the global star formation efficiency, and that the corresponding enriched gas from groups has been 
ejected into intergalactic space, presumably ending up in the form discussed in § §|2.4.1 and p. 4. 2 



2.5. Summary of the Low Redshift Budget 

Our estimates of the baryon budget at low redshift, expressed as fractions of the critical 
Einstein-de Sitter value, are listed in lines 1 to 7 in Table 3. We rewrite the Hubble constant from 
equation (|l]) as H$ = 70hjo km s _1 Mpc -1 to center it on a more likely value. The central values in 
the table are computed from the central values of the parameters entering the estimate. An upper 
bound is obtained by systematically choosing each parameter at the end of its reasonable-looking 
range that maximizes the contribution to the mean mass density, and likewise for the lower bound. 

Line 7a is an estimate of the X-ray-detected plasma around groups (eq. |2^]), line 7b is an 
estimate of the cooler component detected by resonance Lyman a absorption in the trace HI 
fraction (eq. f3lj), and the sum is an estimate of the plasma associated with groups of galaxies. 
Both entries are quite uncertain, and may well have missed a substantial amount of cool plasma. 
Consistent with this, the entry in line 7' from the scaling from the plasma mass in clusters of 
galaxies (eq. |32|| ) is larger. We suspect line 7' is more reliable than the sum of lines 7a and 7b. In 
the sums over the budget we use line 7' for the maximum and central values, and lines 7a and 7b 
instead of 7' for the minimum. 
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The true values for each line are not likely to be close to the maximum in all entries, or close 
to the minima in all entries, but we retain the whole range in quoting the sums in each column in 
line 8 because the sums are dominated by a few entries. 

At hjo = 1 our estimate of the baryons in stars and remnants is £2 = 0.0035 ^0 0014 times the 
critical density. This star mass is only about 17% of our sum over all detected baryons. In our 
budget, in the 2.2// K band, 26% of the light is from disk and irregular galaxies, in agreement 
with the fraction of the star mass. (We use (B — K) = 4.27 for ellipticals and bulges, and 2.78 
for disks and irregulars). In the B band, 61% of the mean luminosity density is from disks and 
irregular galaxies. This component represents only about 25 % of the star mass and only about 
4% of our estimated sum of all observed baryons. 

2.6. Uncounted Components 

We discuss here several mass components that are not counted in our budget because the 
observational evidence about them is still sketchy. 

2.6.1. Warm Plasma in the Voids 

There could be a very large baryon mass in plasma in the voids. Its emission may even have 
been observed: plasma at temperature ~ 2 x 10 6 K fits a component of the diffuse soft X-ray 
background (Wang & McCray 1993). This might be produced locally by a relatively small amount 
of gas in the Galaxy (for a summary of evidence for and limits on hot gas around the Milky Way 
see Moore & Davis 1994), or it could be produced by diffuse extragalactic gas with a mean density 

(fiffli)iGM = OlC^C-^h^ 2 , (33) 

where the emissivity parameter ranges from ( = 1 for solar to £ = 0.1 for primordial abundances, 
and the clumping parameter C = (p? e ) / (n e ) 2 may be significantly greater than unity. Gas at this 
temperature at low redshift has few other detectable effects. For example, the COBE limit on the 
Compton distortion of the microwave background spectrum, y < 1.5 x 10~ 5 (Fixsen et al. 1996), 
is produced by a Hubble length of gas at a mean density Q = (T/10 6 ) - 1 , so the bound on the 
y-distortion is not useful here. On theoretical grounds however it is difficult to see how there could 
be much mass in void plasma. If primeval, it would be surprising that the voids contain little else 
in the form of baryons, especially as detected in diffuse helium absorption at z « 3 (see below). 
If blown in, it would be surprising if the plasma far exceeds the density of stars, and the required 
temperature in this case (to achieve sufficient gas velocity to refill the voids in a Hubble time) 
exceeds 10 7 K, which does produce excessive anisotropy and spectral distortion in the background 
radiation. 
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2.6.2. MACHO s 

The massive compact halo objects (MACHOs) detected as gravitational microlensing events 
could be stellar remnants, that is, baryons. The nature of the MACHOs is not known, and 
there is no secure estimate of their global abundance. However, they do seem to comprise a new 
population not otherwise accounted for. 

Current results from the MACHO collaboration (Alcock et al. 1997) indicate objects with 
mass comparable to that of the Sun (in one model, 0.5^q2-Mq) may account for 20% to 100% of the 
dark mass in a standard spherical halo between the Milky Way and the Large Magellanic Cloud. 
This result must be taken with caution, however, because the experiment measures the MACHO 
mass column in just one direction in one halo and the extrapolation to a global density is subject 
to many uncalibrated assumptions. The Galactic MACHO population may have asymmetries or 
the MACHOs may be concentrated relative to the global dark matter in the halo. For example, 
the actual mass of MACHOs inferred within 50kpc (based on a spherical model, but without 
extrapolating to larger radii) is estimated at 68% confidence to be 13 to 32xlO lo M . This can be 
compared to the mass of the disk (about 6 x 10 10 -M©). A reasonable estimate of the minimum 
global density of MACHOs, viewing them as a new Galactic stellar population, is to assume that 
all disk galaxies have the same ratio of MACHO to disk mass, taking the low end of the estimated 
range; this yields 

^MACHO,min = 2.2J7disks,min = 0.0011/l 7 Q, (34) 

a minor entry in the budget. On the other hand the data are consistent with all of the 
Galactic dark matter being in MACHOs, so a reasonable upper limit derives from assuming that 
MACHOs comprise 100% of the density parameter in gravitation mass in equation (^), giving 
^MACHO.max = 0.25, and making MACHOs the dominant entry in the baryon budget. (At this 
level one expects detectable effects due to microlensing of quasar continuum emission regions, e.g. 
Dalcanton et al. 1994). 

Assessments of this serious uncertainty will be guided by advances in observational constraints 
on the nature and amount of the MACHOs, and perhaps also by advances in understanding 
processes of star formation and death that could produce a substantial mass in baryonic MACHOs 
without leaving an unacceptably large amount of debris. 

2.6.3. Dwarf Galaxies and Low Surface Brightness Galaxies 

The fainter end of the galaxy luminosity function has been a matter of debate for some time, 
but recent work with larger volume surveys (Lin et al. 1996; Zucca et al. 1997) indicates the 
luminosity function is close to flat for at least 5 mag down from L* for the local field galaxies. 
There is some evidence for a rise of the fainter tail right after the shoulder of L* (Ellis et al. 1996), 
but this is visible only in a high z sample selected in the blue band, for which small galaxies 
show extraordinary star formation activity without contributing much to the baryon budget. This 
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means the contribution of dwarf galaxies to the baryon budget is small, and the estimate from 
the local luminosity density likely is more secure. A rise in the luminosity function somewhat 
below L* for cluster members also has been reported. (For the most recent literature see Phillipps 
et al. 1997). The contribution to the mass density from these cluster dwarf galaxies cannot be 
dominant unless this sharp rise continues to very low mass, which is not likely since it would 
violate measurements of cluster surface brightnesses. 

The review by Bothun, Impey, & McGaugh (1997) indicates the importance of galaxies 
with surface brightnesses lower than the more readily detectable "normal" galaxies. It is not 
clear however that low surface brightness (LSB) galaxies make a significant contribution to the 
baryon budget. LSB galaxies contribute to the budget by their stars and star remnants, neutral 
gas, and plasma. The 21-cm observations by Briggs (1997) indicate there is not a significant 
contribution to the mean density of atomic hydrogen from gas-rich LSB galaxies in the field at 
distances <^ 10ft, -1 Mpc. A direct constraint on diffuse plasma around LSB galaxies is much 
more difficult, but since LSB galaxies avoid the voids defined by the more visible galaxies this is 
the problem of counting the plasma concentrated around galaxies outside the great clusters, as 
summarized in lines 7a, 7b, and 7'. As we have noted, diffuse baryons initially associated with 
LSB galaxies that end up as cluster members would be counted in the X-ray measurements of 



(^Hii)ci (eq. pM ) and in the scaled value for the plasma in the field (eq. [32] and line 7' in Table 3). 
The mean luminosity density in the LSB sample considered by Sprayberry et al. (1997) is 15% 
of our adopted value (eq. The mean luminosity density from all stars not shrouded by dust, 
including those in systems with luminosities or surface brightnesses below detection thresholds, 
is constrained by measurements of the extragalactic contribution to the sky surface brightness. 
Absolute measurements in progress (Bernstein 1997) are capable of reaching ~ 3C (eq. [Q]), but 
differential measurements of background fluctuations already constrain plausible new populations 
to contribute much less than £ (Dalcanton et al. 1997, Vogeley 1997). We cannot exclude the 
possibility that there is a significant mass in brown dwarfs or baryonic MACHOs in LSB galaxies, 
but we can note that since many of these galaxies seem to be in early states of evolution they 
would not seem to have had much opportunity to have sequestered mass in dark stars, and that 
the integrated baryon content in present-day LSB galaxies likely is counted in the measures of 
diffuse gas at redshift z ~ 3, as discussed next. 



3. The Baryon Budget at z ks 3 

There has been significant progress in understanding stellar populations at high redshift (e.g., 
Madau et al. 1996), but many issues remain open. It is fortunate for our purpose that stars are 
subdominant in the budget at the present epoch and likely are even less important at high redshift, 
so we can concentrate on the constraints on diffuse gas from quasar absorption line spectra. 

The neutral hydrogen at z ~ 3 is predominantly in the high column density damped Lyman-a 
absorbers (DLAs; Lanzetta, Wolfe, & Turnshek 1995 and references therein). The amount of 
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neutral hydrogen in the DLAs increases with increasing redshift back to z ~ 2. Our adopted 
value for the the density parameter in this form at 2 < z < 3 is from Storrie-Lombardi, Irwin, & 
McMahon (1996): 

^atomic(^ = 3) = (0.0013 ± 0.0003, 0.0020 ± 0.0007)(/i/70) _1 . (35) 

The two estimates are for $7 = and $7 = 1 respectively, and include the mass of the accompanying 
helium. Because these systems are optically thick to ionizing radiation there is no correction for 
an ionized fraction in the neutral gas. (The significant amount of mass in plasma in HII regions 
around young stars or in regions exposed to the intergalactic ionizing radiation are included in the 
forest component below). Equation (|35|) could be an underestimate if extinction by dust in the 
gas suppressed the selection of quasars behind high column density absorbers (Fall and Pei 1993) 
or an overestimate if gravitational lensing enhanced selection of lines of sight through dust-free, 
gas-rich absorbers. 

Yet another uncertainty is the residency time of gas in DLAs. The large velocities in the 
DLAs (Prochaska & Wolfe 1997) indicate the cloud masses could be depleted by winds only if 
there were considerable energy input from supernovae. The more likely scenario is that the density 
parameter in DLAs is decreasing at z <^ 2 because the HI is being converted into stars. The HI 
mass in DLAs at z = 3 is about half that in present-day stars (lines 1, 2, 3, and 9 in Table 3). 
That could mean there is a significant mass in stars in DLAs and other young galaxies at z = 3 
and/or that intergalactic matter still is settling onto protogalaxies at z = 3. However, these issues 
do not affect the budget regarded as a snapshot of conditions at z = 3. 

The dominant baryonic mass component is the Lyman-a forest gas, detected by the trace 
neutral hydrogen in plasma that fills space as a froth at z = 3. The density can be estimated 
from CDM model simulations that give good fits to the observations taking into account the 
distributions of cloud shapes, sizes, flow velocities, and temperatures, although these estimates 
are still sensitive to the uncertain flux of ionizing radiation. Rauch et al. (1997) conclude that the 
baryon density parameter needed to correctly reproduce the statistical absorption (mostly due to 
clouds at HI surface density ~ 10 13±1 cm -2 ) is $7hii/i 2 > 0.017 — 0.021, depending on cosmological 
model. Weinberg et al. (1997) quote a lower limit of 0.0125/i~ 2 . Zhang et al. (1997), including a 
self-consistent analysis of the ionizing spectrum, arrive at a range for h = 0.5 of 0.03 < $7^ < 0.08, 
with about half of this in the forest clouds. Smaller densities may be possible however because the 
same absorption can be produced by a higher HI fraction, caused by higher density contrast and 
lower gas entropy (Wadsley and Bond 1996, Bond and Wadsley 1997). Because of the unresolved 
issues we assign the current estimates a "B" grade. We adopt in Table |3| the range given by Zhang 
et al. (1997), scaling by /i~ 3 / 2 . For the lower limit we include not the total density but just that in 
the forest clouds; this is appropriate for consistency, since the simulations predict a larger fraction 
of cooled baryonic matter than we infer from observations. 

Cool plasma between the forest clouds has lower density and hence a lower neutral fraction 
and very low HI Lyman-a resonance optical depth. The amount of plasma in this form is 
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best probed by measurements of resonance absorption by the most abundant absorbing ion, 
singly ionized helium (Jakobsen et al. 1994; Davidsen, Kriss, & Zheng 1996.) High resolution 
HST/GHRS quasar absorption line spectra of two quasars now permit the separation of the more 
diffuse component of the Hell resonance absorption from the component in the HI Lyman-a forest 
clouds (Hogan, Anderson, & Rugers 1997, Reimers et al. 1997). The upper bound on intercloud 
gas density is derived based on the maximum permitted mean Hell Lyman-a optical depth allowed 
after subtracting the minimal contribution from the detected HI Lyman-a forest clouds, while 
adopting the hardest ionizing spectrum allowed by the data, with opposite assumptions leading to 
the lower bound. (Note that the upper bound is on photoionized cool gas; hot gas with thermally 
ionized helium is not constrained by absorption, but is even less plausible at z « 3 than at z ~ 0.) 
These limits are consistent with predictions from numerical simulations of CDM models, in which 
most of the gas is clumped in redshift space by this time (Croft et al. 1997, Zhang et al. 1997, Bi 
and Davidsen 1997.) Although in principle the helium bounds are sound the present results are 
given a "B" grade because they are as yet based on only two quasars. 

4. Nucleosynthesis 

The observationally successful theory of the origin of the light elements by nucleosynthesis at 
redshift z ~ 10 9 predicts the mean baryon density in terms of the primeval element abundances 
(as reviewed extensively in the literature; for example Walker et al. . 1991, Copi et al. 1995, Hata 
et al. 1997). We consider here standard homogeneous nucleosynthesis predictions for abundances 
as a function of the baryon-to-photon ratio r\ = 10~ 10 ?7io, where the present baryon density is 
^baryon ^70 = 7.45 x 10 — 3 77io- In this context the strongest constraints on fibaryon derive from the 
abundances of helium and deuterium. 

The primordial deuterium abundance is still uncertain but we can quote reliable upper and 
lower bounds. A conservative lower bound, (D/H) p > 2 x 10~ 5 , comes from many sources — the 
Jovian atmosphere (e.g. Niemann 1996), the interstellar medium (Ferlet & Lemoine 1996, Linsky 
et al. 1995), and quasar absorption lines (Tytler, Fan & Buries 1996), together with the fact that 
no source outside the Big Bang is known to produce deuterium significantly. An upper bound, 
(D/H) p < 2 x 10~ 4 , comes from several measurements in metal-poor quasar absorption systems 
(Songaila et al. 1994; Carswell et al. 1994; Webb et al. 1997; Songalia, Wampler, & Cowie, 
1997). Although the identification as deuterium in these systems is disputed (Tytler, Buries, &; 
Kirkman 1997), the lack of higher detected values, and the fact that significant D destruction 
would normally be expected to produce significant metal enrichment, makes this a robust upper 
limit. These yield the limits 1.7 < 7710 < 7.2. The lower bound is used for the minimum value of 
^baryon i n hne 14 in Table ||| For comparison we include the central value favored by Buries and 
Tytler (1997), (D/H) p = 3.4 x 1(T 5 , which yields r?i = 5.1 or ft baryon = 0.039/ifo 2 . 

The helium abundance Y is well measured in nearby galaxies (e.g., Pagel et al. 1992; 
Skillmann et al. 1994; Izotov, Thuan, & Lipovetsky 1997). The primordial value Y p derived 
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from these observations depends on models of stellar enrichment, but the present datasets yield 
a nearly model-independent 2a Bayesian upper limit Y p < 0.243 (Hogan, Olive, & Scully 1997). 
This corresponds to t]iq = 3.6, which we use for the maximum value in line 13 in Table ||; it is 
about half the upper limit from deuterium. Most current studies (e.g. Olive & Steigman 1995) are 
consistent with the central value Y p ~ 0.23 used in the table. 



5. Discussion 



After a brief comparison with earlier work (§ [54]) we offer in § 5.2 an assessment of the major 
uncertainties in our budget at low redshift, with emphasis on the more indirect arguments we hope 
constrain the budget of baryons in forms that are unobservable in practice. We review the test 



from the theory for the origin of the light elements in § [h3, and in § 5^ we compare these results 
with measurements of the diffuse baryons at z ~ 3. Some concluding remarks are offered in § 5.5 



5.1. Comparison with Earlier Estimates 

Let us compare our estimates with those by Persic & Salucci (1992), Gnedin & Ostriker 
(1992), and Bristow & Phillipps (1994). These three groups considered baryons in stars (often 
classified into those in early and late types of galaxies) and in hot X-ray emitting gas around 
clusters. The baryon abundance in stars estimated by Gnedin & Ostriker is in good agreement 
with ours, and that by Persic & Salucci is close to our minimum estimate. On the other hand, the 
estimate of Bristow & Phillipps is 6 times ours. The same pattern applies to baryons in X-ray 
emitting gas: Persic &: Salucci obtained 1/3 our central value and Bristow &; Phillipps found twice 
our estimate. As a result Persic &: Salucci find that the total amount of baryons is substantially 
less than that expected from nucleosynthesis, while Bristow &: Phillipps find the nucleosynthesis 
value is saturated by stars plus hot X-ray emitting gas alone. Our attempt to take account of a 
broader range of constraints and states of baryons, including plasma in groups inferred but not 
directly observed, has led us to conclude that, apart from baryons that may have been sequestered 
at very high redshift, the bulk are either directly visible or in forms which are well constrained 
and accounted for with plausible extrapolations. 



5.2. Observable and Unobservable Baryons 

A baryon budget must be informed by an understanding of what is reasonable and sensible 
within standard models for the physics, astronomy, and cosmology. We argue here that these 
considerations coupled to the observations available now or within reach significantly limit the 
possibilities for large gaps in the budget. 
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Chemically Bound Baryons. Baryons confined by chemical bonds, as in comets, are dark 
(apart from evaporation) at any significant distance, but we know their contribution to the baryon 
mass within the luminous parts of normal galaxies is not large: the mass in stars accounts for 
the dynamics. In the standard cosmology only a small fraction of primordial hydrogen converts 
to molecular form. Some scenarios envisage a significant mass in molecular clouds (Combes and 
Pfenniger 1998), but it seems more likely to us that the bulk of this material would by now have 
formed into stars, as it has in the Galaxy. Baryons could be chemically bound to heavier elements, 
but it is hard to see how this could be a significant mass component because it would require 
prodigious heavy element production in closely confined and closed systems and the radiated 
nuclear binding energy would appear in the background radiation. Besides these arguments we 
have the direct limits on the density of the plentiful but inert gas helium which would separate 
from any condensed component. In short, if the existence of a significant mass in intergalactic 
comets could be established it would contradict the standard cosmological model. The far more 
credible proposition is that there is not a significant mass fraction in chemically bound baryons. 

Dark Galaxies. The starlight from baryons gravitationally bound in known forms of low 
surface brightness (LSB) galaxies is only about 15% of our adopted total. Thus if normal galaxies 
were an adequate guide to their compositions the known LSB galaxies would be insignificant 
reservoirs of baryons in stars. Discovery of reason to suspect the contrary would be of great 
interest, of course. The starlight in as yet undetected still lower surface brightness galaxies is 
constrained by the absolute surface brightness of the sky and fluctuations thereof; a significant 
contribution already appears unlikely from the smoothness of the background in deep exposures. If 
measurements indicated the mean optical extragalactic surface brightness is larger than expected 
from known galaxies it certainly would give reason to think we have missed an important baryon 
component. 

If the surface brightness showed no such anomaly one could still imagine there are dark 
galaxies with star populations dominated by brown dwarfs, but there are several arguments against 
the idea. There are high success rates in optical identification of systems of normal galaxies with 
the objects responsible for gravitational lensing events, for quasar absorption line systems, and for 
X-ray sources. Even if massive dark galaxies were assumed not to be capable of significant lensing 
they could contain diffuse baryonic halos, so in this scenario we would need to explain why they do 
not cause a significant rate of identification failures for quasar absorption lines and X-ray sources. 
One way assumes the dark galaxies avoid the concentrations of normal galaxies and the associated 
plasma, so they do not acquire plasma halos, but there are two counter arguments. First, in this 
picture there surely would be intermediate cases, dim massive galaxies that avoid the normal ones, 
and there is no evidence of them. Second, the familiar morphology-density relation in observed 
galaxies goes the other way: gas-rich galaxies prefer lower density environments, gas-poor anemic 
spirals denser regions. 

Dust-Shrouded Stars. Populations of stars shrouded in dust are constrained by the integrated 
far infrared background of reradiated starlight, which is observed to have an integrated flux 
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about twice that from optical galaxies as estimated from the Hubble Deep Field and is thus the 
repository of most of the nuclear energy released by stars (Schlegel, Finkbeiner and Davis 1998; 
Hauser et al. 1998). The total mass in baryons in stars even allowing for a maximal shrouded 
population must however still be small. It is also interesting to compare the estimated total energy 
density with the corresponding global production of heavy elements. Taking the total bolometric 
intensity of starlight to be about 50 nW m~ 2 ster _1 (based on Schlegel et al.'s fitted value 32 ± 13 
at 140/xm), the energy density is 

u = 2 x 1(T 14 (1 + z)erg cm" 3 = 1 x 10~ 8 (1 + z)MeV cm" 3 (36) 

where the redshift factor corrects for energy lost after the mean epoch of emission. Each 
nucleon that is burned to helium releases 25 MeV in heat, rising to a total of about 30 MeV 
per nucleon converted from hydrogen to heavy elements. On dividing u by 30 MeV we get the 
mean nucleon number density in heavy elements produced in the production of the background 
light, rihcavy ~ 4 x 10~ 10 (1 + z)cm~ 3 . Our central value for the baryon number density is 
"-baryons = 11 x 10~ 7 cm~ 3 ; the ratio is the mass fraction in heavy elements to make the observed 
background, 

Z ~ 0.004(1 +z). (37) 

If the bulk of the radiation were produced at z = 2 it would mean Z is about 1 percent. Though 
not a precise constraint, it is a significant check that the number roughly agrees with the observed 
metallicity of baryonic material. 

Black Holes. Baryons sequestered in black holes in the luminous parts of normal galaxies 
are known to be subdominant to the baryons in stars. It is quite unreasonable to imagine whole 
galaxies have been lost to relativistic collapse. Gas clouds at the Jeans length at decoupling may 
be susceptible to relativistic collapse but only a tiny fraction of material has small enough angular 
momentum to form black holes, and even then it is hard to imagine an appreciable fraction of the 
cloud mass could collapse before the remainder is blown apart by radiation from the collapsing 
fraction. Far before recombination even a tiny collapsing fraction yields a significant mass fraction 
today, although it requires either extremely large amplitude perturbations in 77, which tends to 
adversely affect the model for light element production (e.g., Kurki-Suonio, Jedamzik, & Mathews 
1996), or extreme fine tuning of the spectrum and amplitude of small-scale adiabatic perturbations 
(Carr 1994). In these scenarios baryons may be sequestered in compact objects at redshifts 
z ^ 10 10 that are in effect nonbaryonic for the purpose of this paper. 

Baryonic MACHOs. Debris is a key issue for MACHOs interpreted as star remnants. In 
known processes of formation of white dwarfs and neutron stars, winds and explosions disperse 
most of the original star mass in diffuse debris. Thus a 2M Q star leaves a 0.5M Q remnant and 
1.5Mq in debris, and the ejected fraction is larger for larger stars. The lensing observations are 
consistent with a density parameter Q ~ 0.25 in MACHOs, which would be a very considerable 
entry in the budget, but we would have to explain what happened to the debris. The debris might 
be recycled and a large fraction ultimately sequestered in many generations of MACHOs, but 
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recycling seems unlikely in the low density halos of galaxies. Debris still present as diffuse matter, 
with density parameter f2 > 0.014, would violate constraints previously discussed on diffuse matter 
in groups and most models would also generate excessive metals. The most plausible form of 
baryons is brown dwarfs, though this would require a second peak in the global IMF. Brown dwarf 
masses are below the estimated mass range of the observed MACHOs, but this could be the result 
of an unfortunate distribution of velocities of the MACHOs in the direction of the LMC. Absent 
the discovery of observational evidence for the brown dwarf picture, or the demonstration that 
the formation of star remnants need not disperse much debris in observable forms, we suspect the 
MACHO contribution to the baryon budget is subdominant to the mass in plasma around groups. 

Plasma Around Field Galaxies. There is a significant baryon density in stars in normal 
galaxies, and a still larger density in plasma around the galaxies. The latter is the largest and most 
uncertain entry in our baryon budget at low redshift in Table [| The estimate in line (7') assumes 
the spheroid star-to-gravitational mass ratio is the same in clusters and the field. If we have 
missed a significant mass in spheroids, as in LSB galaxies, and the missing fraction is the same 
in clusters and the field, it does not affect the estimate. If spheroid star production were more 
efficient in clusters it would mean that the estimated group mass is too small and hence line (7') 
is too small. If plasma were ejected (more than spheroid stars) from clusters during assembly, 
line (7') would be an overestimate. But neither of these effects could be large without upsetting 
the condordance between equations (24), (25) and (27). Comparison of the two estimates, line (7') 
and the sum of lines (7a) and (7b), shows a factor of two difference, but it is easy to imagine 
that this is because the X-ray observations still miss considerable warm plasma and we have not 
yet adequately modeled the cooler low surface density clouds detected by absorption lines at low 
redshift. 

Baryons in the Voids. There are galaxies in the low density regions, or voids, defined by 
normal galaxies. Galaxies in low density regions tend to be the later Hubble types, and so are more 
readily observable than their counterparts in concentrations. Gas enrichment aside, no known 
type of object prefers the voids. This includes dwarf and irregular galaxies, low surface brightness 
galaxies, and the plasma clouds detected from absorption lines. The reasonable presumption is 
that plasma clusters the way all other observed baryons do, meaning there is not much mass in 
plasma in the voids. The possible exception is gas with high enough pressure to resist gravitational 
draining of the voids, but if the density of such void plasma were comparable to that estimated in 
groups it would have appeared already as excessive helium absorption at z ~ 3. 

The lack of reasonable-looking alternatives leads us to conclude that a fair accounting of the 
the baryons is possible because most are in states that can be observed or reliably constrained by 
more indirect arguments. We now argue that there is a reasonable case for the net baryon density 
parameter 

^baryon = 0.021 ± 0.007, (38) 

close to the central value of the sum in line 8 in Table ||. 
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5.3. The Nucleosynthesis Check 

The central density in equation (^) with hjo = 1 corresponds to 7710 = 2.8. At this value 
of the baryon-to-photon ratio the standard model for light element nucleosynthesis predicts 
yj,( 4 He) ~ 0.238 and D/H~ 9 x 10~ 5 . This is consistent with observations of the helium 
abundance, if somewhat above the central value, and with current observations of the deuterium 
abundance, if somewhat toward the high end of the range of estimates (Songaila 1997). The mean 
of the upper limit from helium and the lower limit from deuterium is 7710 = 2.7 which we adopt as 
the central value in line 14 of Table 3. 

If the baryon density were at the low end of our range of estimates in line 8 in Table 3, 
^baryon = 0.007 (t]iq = 1.0 at h?Q = 1), the standard model prediction for the deuterium abundance 
would be unacceptably large. This could be remedied by assuming a significant baryon mass has 
been sequestered in the MACHOs, but for the reasons explained above we are inclined to suspect 
rather that the adopted minima in lines 7a and 7b in Table 3 seriously underestimate the mass 
in plasma. The baryon density at the high end of our estimates is ^baryon = 0.041 (t]iq = 5.5 at 
^•70 = 1)- This density is consistent with the bound from the deuterium abundance but it predicts 
Y p = 0.247, well outside the accepted bound on the primeval helium mass fraction. In our opinion 
the likely interpretation is that our upper limit in line 8 of Table 3 also is overly conservative. 

Our interpretation of the baryon budget assumes the standard homogeneous model for 
primordial nucleosynthesis. An inhomogeneous primeval entropy per baryon, with the appropriate 
coherence length (Jedamzik &: Fuller 1995), would require serious attention if further observations 
confirmed that there are significant variations in the deuterium abundances in high redshift HI 
clouds or that the helium abundances in some dwarf galaxies are lower than that predicted by the 
global baryon budget applied to the homogeneous nucleosynthesis model. 

5.4. The Check from the Budget at z f» 3 

Stars are subdominant in the budget at low redshift and likely are even less significant at 
redshift z ~ 3, so it makes sense to compare the sum of the high redshift diffuse gas components 
in Table 3 to the sum at low redshift and the prediction of light element nucleosynthesis. 

The sum of the central values for diffuse baryons at z ~ 3 is Q = 0.04. This is consistent with 
the maximum sum in our low redshift budget but again contradicts the standard model for helium 
abundance. The sum of the minimum values at z ~ 3 is within the error flags in equation ( |38|) and 
the standard model for the light elements. Thus there is no contradiction with equation (|3S|), but 
advances in the observations of the intergalactic medium at high redshift may provide a critical 
test. 
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5.5. Concluding Remarks 

We have emphasized three themes. First, a robust estimate of the present baryon budget is of 
central importance in serving to clarify and sharpen issues of research in present-day cosmology. 
Second, there is no guarantee that physical processes in the formation and evolution of structure 
will have conspired to place a significant fraction of the baryons in forms we can hope to detect 
by their emission or absorption of radiation or otherwise infer in a convincing way. Third, the 
developing evidence is that Nature in fact has been kind: we may actually be able to arrive at a 
close to complete budget of the baryons in the observable universe. 



The value for Ob ar yon i n equation (38) suggests or is consistent with three lines of ideas. 
First, it suggests baryonic MACHOs make an insignificant contribution to the global budget. 
Perhaps the MACHOs are non-baryonic entities formed in the early universe, or perhaps like the 
stellar components of the Milky Way galaxy they are more concentrated to the inner halo than is 
nonbaryonic dark matter. The second idea is that there are relatively few baryons in the voids 
defined by the galaxies, likely because the baryons were swept out of the voids by gravity. Gravity 
would gather low pressure dark matter with the baryons and galaxies, implying that the mass 
density in matter than can cluster is well below the Einstein-de Sitter value. This is consistent 
with our estimate of the total density parameter in equation (p4|). The third idea is that the bulk 
of the forest baryons at z = 3 have ended up now in warm plasma around the galaxies outside the 
great clusters. The mass in neutral hydrogen in the DLAs is about half that needed to account 
for the present-day mass in stars. We know some mass already is in stars at z = 3; formation of 
the rest of the star mass at z < 3 would be a modest drain on the mass in the forest. 

An interesting feature of the budget is that all stars and their remnants, together with cold 
phases of collapsed gas, comprise a small fraction even of the small fraction of the matter which 
is baryonic (Gnedin & Ostriker 1992; Persic & Salucci 1992). Just the 0.3% of the baryons in 
irregular galaxies shine brightly enough at z = 0.4 to 1 to make an important contribution to the 
rapid increase of B band number counts of galaxies with increasing apparent magnitude. It would 
not be difficult in these circumstances to imagine that the distribution of the starlight has little to 
do with the mass. But the consistency of the estimates of dynamical estimates of the total density 
parameter, as in equations (p4) , ( pq) , (p?|) suggests starlight nevertheless is a useful tracer of mass 
on the scale of the distance between galaxies. 

The baryon budget still is quite uncertain: the sums in line 8 in Table || for the maximum 
and minimum low redshift budget differ by a factor of six. We suspect this is largely a result of 
an underestimate in line (7a) for plasma in groups; improved information on soft X-ray emission 
would be of considerable help. Perhaps improvements in the observational constraints will 
confirm the tentative concordance in equation (^); perhaps better observations will reveal an 
inconsistency that shows we have to reconsider some aspect of the standard concepts of astronomy 
and cosmology. 



-24- 



We are grateful to J. Bahcall for detailed discussions of the local star luminosity function 
and to R. Carlbcrg, J. Stadel, C. Stubbs, and E. Vishniac for useful comments. This work was 
supported in part by a Grant-in- Aid of the Ministry of Education of Japan (09640336) and 
the Japan Society for the Promotion of Science, by NASA and the NSF at the University of 
Washington, and by the NSF at Princeton University. M.F. acknowledges support from the Fuji 
Xerox Corporation at Princeton. 



-25- 



REFERENCES 

Alcock, C. et al. 1997, ApJ, 486, 697 

Bahcall, J. N. & Soneira, R. M. 1980, ApJS, 44, 73 

Bahcall, N. A. & Cen, R. 1993, ApJ, 407, L49 

Bahcall, N. A., Fan, X., & Cen, R. 1997, ApJ Letters, in press 

Bahcall, N. A., Lubin, L., & Dorman, V. 1995, ApJ, 447, L81 

Bernstein, R. A. 1997, Ph. D. thesis, Caltech 

Bi, H. & Davidsen, A. F. 1997, ApJ 479, 523 

Binney, J. & Tremaine, S. 1987, Galactic Dynamics (Princeton: Princeton University Press) 

Bond, J. R. &: Wadsley, J. W. 1997, To appear in Structure and Evolution of the IGM from QSO 
Absorption Line Systems, Proc. 13th IAP Colloquium, ed. P. Petitjean and S. Chariot, 
Nouvelles Frontieres, Paris 

Bothun, C, Impey, C, & McGaugh, S. 1997, PASP, 109, 745 

Briggs, F. H. 1997, ApJ, 484, L29 

Bristow, P. D. & Phillipps, S. 1994, MNRAS, 267, 13 

Buries, S. & Tytler, D. 1997, preprint 

Carlberg, R. C, Yee, H. K. C, Ellingson, E., Abraham, R., Gravel, P., Morris, S., & Pritchet, C. 
J. 1996, ApJ, 462, 32 

Carlberg, R. G. et al. 1997, preprint 

Carr, B. J. 1994, ARA& A, 32, 531 

CarsweU, R. F., Weymann, R. J., Cooke, A. J., k Webb, J. K. 1994, MNRAS, 268, LI 
Chariot, S., Worthey, G., & Bressan, A. 1996, ApJ, 457, 625 
Combes, F. and Pfenniger, D. 1998, preprint 

Copi, C. J., Schramm, D. N., & Turner, M. S. 1995, Science, 267, 192 

Croft, R. A. C, Weinberg, D. H., Katz, N., & Hernquist, L. 1997, ApJ, in press 

da Costa, L. N. et al. 1994, ApJ, 424, LI 

Dalcanton, J. J., Canizares, C, Granados, A., Steidel, C. C, & Stocke, J. T. 1994, ApJ 424, 550 

Dalcanton, J. J., Spergel, D. N., Gunn, J. E., Schmidt, M., k, Schneider, D. P. 1997, AJ, in press 
( |a^tro-ph/9705088D 

Davidsen, A. F., Kriss, G. A., & Zheng, W. 1996, Nature, 380, 47 

de Vaucouleurs, G. & Pence, W. D. 1978, AJ, 83, 1163 

Dressier, A. 1980, ApJ, 236, 351 



-26- 



Efstathiou, G., Ellis, R. S., & Peterson, B. A. 1988, MNRAS, 232, 431 
Ellis, R. S. 1997, ARA& A, to be published 
Ellis, R. S. et al. 1996, MNRAS, 280, 235 

Fabricant, D., Beers, T. C, Geller, M. J., Gorenstein, P., Huchra, J. P., & Kurtz, M. J. 1986, ApJ, 
308, 530 

Fall, S. M. & Pei, Y. C. 1993, ApJ, 402, 479 

Ferlet, R. and Lemoine, M. 1996, ASP Conference Series, 99, 78 

Fixsen, D. J., Cheng, E. S., Gales, J. M., Mather, J. C, Shafer, R. A., & Wright, E. L. 1996, ApJ, 
473, 576 

Fukugita, M., Hogan, C. J., & Peebles, P. J. E. 1996, Nature, 381, 489 

Fukugita, M., Shimasaku, K., & Ichikawa, T. 1995, PASP, 107, 945 

Garnavich, P. M. et al. 1998, ApJ, in press 

Gnedin, N. Yu & Ostriker, J. P. 1992, ApJ, 400, 1 

Gould, A., Bahcall, J. N., & Flynn, C. 1996 ApJ, 465, 759 

Hata, N., Steigman, G., Bludman, S., & Langacker, P. 1997, Phys Rev, 55, 540 

Hauser, M. et al. 1998, preprint 

Hogan, C. J., Anderson, S. F., & Rugers, M. H. 1997, AJ, 113, 1495 
Hogan, C. J., Olive, K. A., & Scully, S. T. 1997, ApJL, in press 
Hughes, J. P. 1989, ApJ, 337, 21 

Izotov, Y. I., Thuan, T. X., & Lipovetsky, V. A. 1997, ApJS, 108, 1 

Jakobsen, P., Boksenberg, A., Deharreng, J. M., Greenfield, P., Jedrzejewski, R., & Paresce, F. 
1994, Nature 370, 35 

Jedamzik, K. & Fuller, G. M. 1995, ApJ, 452, 33 

Kent, S. M. 1985, ApJS, 59, 115 

Kormendy. J. 1986, in Structure and Dynamics of Elliptical Galaxies (IAU Symposium 127), ed. 
by T. de Zeeuw (Dordrecht: Reidel), p. 17 

Kurki-Suonio, H., Jedamzik, K., & Mathews, G. J. 1996, preprint 

Lanzetta, K. M., Wolfe, A. M. & Turnshek, D. A. 1995, ApJ, 440, 435 

Lilly, S. J., Tresse, L., Hammer, F., Crampton, D., & Le Fevre, O. 1995, ApJ, 455, 108 

Lin, H., Kirshner, R. P., Schechter, S. A., Landy, S. D., Oemler, A., Tucker, D. L. & Schechter, P. 
L. 1996, ApJ, 464, 60 

Linsky, J. L., Diplas, A., Wood, B. E., Brown, A., Ayres, T. R., & Savage, B. D. 1995, ApJ 451, 
335 



-27- 



Loveday, J., Peterson, B. A., Efstathiou, G., k Maddox, S. J. 1992, ApJ, 390, 338 

Madau, P., Ferguson, H. C, Dickinson, M. E., Giavalisco, M., Steidel, C. C, k Fruchter, A. 1996, 
MNRAS, 283, 1388 

Marzke, R. O., Huchra, J. P., k Geller, M. J. 1994, ApJ, 428, 43 

Meisels, A. k Ostriker, J. P. 1984, AJ, 89,1451 

Mihalas, D. k Binney, J. 1981, Galactic Astronomy (San Francisco: Freeman) 
Moore, B. k Davis, M. 1994, MNRAS, 270, 209 
Morgan, W. W. 1962, ApJ, 135, 1 

Mulchaey. J. S., Davis, D. S., Mushotzky, R. F., k Burstein, D. 1996, ApJ, 456, 80 

Myers, S. T., Baker, J. E., Readhead, A. C. S., Leitch, E. M., k Herbig, T. 1997, ApJ 485, 1 

Niemann, H. B. 1996, Science, 272, 846 

Olive, K. A. k Steigman, G. 1995, ApJS, 97, 49 

Pagel, B. E. J., Simonson, E. A., Terlevich, R. A., k Edmunds, M. 1992, MNRAS, 255, 325 
Peacock, J. 1997, MNRAS, 284, 885P 
Peebles, P. J. E. 1986, Nature, 321, 27 
Perlmutter, S. et al. 1998, Nature 391, 51 
Persic, M. k Salucci, P. 1992, MNRAS, 258, 14p 

Phillipps, S., Parker, Q. A., Schwartzenberg, J. M. k Jones, J. B. 1997, ApJ, L59 
Postman, M. k Geller, M. J. 1985, ApJ, 281, 95 
Prochaska, J. X. k Wolfe, A. M. 1997, ApJ, in press 
Rao, S. k Briggs, F. 1993, ApJ, 419, 515 

Rauch, M., Miralda-Escude, J., Sargent, W. L. W., Barlow, T. A., Weinberg, D. H., Hernquist, L, 
Katz, N., Cen, R., k Ostriker, J. P. 1997, ApJ in press 

Reimers, D., Koehler, S., Wisotzki. L., Groote, D., Rodriguez-Pascual, P., k Wamsteker, W. 
1997, AA 327, 890 

Renzini, A. 1997, ApJ, 488, 35 

Roberts, M. S. k Haynes, M. P. 1994, ARA& A, 32, 115 
Schechter, P. L. k Dressier, A. 1987, AJ 94, 563 
Shimasaku, K. k Fukugita, M. 1997, ApJ in press 
Shull, M., Stocke, J. T., and Penton, S. 1996, AJ, 111, 72 

Skillman, E., D., Terlevich, R. J., Kennicutt, R. C, Garnett, D. R., k Terlevich, E. 1994, ApJ 
431, 172 



-28- 



Schlegel, D. J., Finkbeiner, D. P. k Davis, M. 1998, preprint 
Songaila, A. 1997, preprint 

Songaila, A., Cowie, L. L., Hogan, C, k Rugers, M. 1994, Nature, 368, 599 

Songaila, A., Wampler, E. J., k Cowie, L. L. 1997, Nature, 385, 137 

Sprayberry, D., Impey, C. D., Irwin, M. J., k Bothun, G. D. 1997, ApJ, 482,104 

Storrie-Lombardi, L J., Irwin, M. J., k McMahon, R. G. 1996, MNRAS, 283, L79 

Tinsley, B. M. 1980, ApJ, 241, 41 

Tinsley, B. M. 1981, MNRAS, 194, 63 

Tinsley, B. M. k Gunn, J. E. 1976, ApJ, 203, 52 

Tytler, D., Buries, S,. k Kirkman, D. 1997, ApJ, in press 

Tytler, D., Fan, X. M., k Buries, S. 1996, Nature, 381, 207 

van der Marel, R. 1991, MNRAS, 253, 710 

Vogeley, M. S. 1997, ApJ, submitted 

Wadsley, J. k Bond, J. R. 1996, to appear in "Computational Astrophysics", Proc. 12th Kingston 
Conference, Halifax, Oct. 1996, ed. D. Clarke k M. West (PASP) 

Walker, T. P., Steigman, C, Schramm, D. N., Olive, K. A., k Kang, H. 1991, ApJ, 376, 51 

Wang, Q. D. k McCray, R. 1993, ApJ, 409, L37 

Webb, J. K., Carswell, R. F., Lanzetta, K. M., Ferlett, R., Lemoine, M., Vidal-Madjar, A., k 
Bowen, D. V. 1997, Nature 388, 250 

Weinberg, D., Miralda-Escude, J., Hernquist, L., k Katz, N. 1997, ApJ, submitted 

White, D. k Fabian, A. C. 1995, MNRAS, 273, 72 

White, S. D. M., Navarro, J. F., Evrard, A. E. k Frenk, C. S. 1993, Nature, 366, 429 

Whitmore, B. C. k Kirshner, R. P. 1981, ApJ, 250, 43 

Young, J. S. k Scoville, N. Z. 1991, ARA& A, 29, 581 

Zhang, Y., Anninos, P., Norman, M. L., k Meiksin, A. 1997 ApJ 485, 496 

Zucca, E. et al. 1997, AA 326, 477 

Zwaan, M., Briggs, F., Sprayberry, D., k Sorar, E. 1998, ApJ, in press 



This preprint was prepared with the AAS IATgX macros v4.0. 



-29- 



Table 1. Galaxy Parameters 





E 


SO 


Sab 


Sbc 


Scd 


Irr 


bulge fraction n(r) 


1 


0.75 


0.40 


0.24 


0.10 





bulge fraction n(B) 


1 


0.64 


0.33 


0.16 


0.061 





morphology fraction /x(B) 


0.11 


0.21 


0.28 


0.29 


0.045 


0.061 


HI typical mass (1O 9 M ) 





0.64 


2.5 


4.0 


2.9 


1.8 


H 2 /HI 




2.5 


1.0 


0.63 


0.25 


0.06 



Table 2. Luminosity Function Parameters 



Author 


M*(B) 


a 




£(B) b 


Efstathiou et al. (1988) 


-19.68 


-1.1 


1.56 


1.93 


Loveday et al. (1992) 


-19.50 


-0.97 


1.40 


1.35 


da Costa et al. (1994) 


-19.5 


-1.20 


1.5 


1.71 


Marzke et al. (1994) 


-18.8 


-1.0 


4.0 


2.1 


Lilly et al. (1995) 


-19.36 


-1.03 


2.48 


2.18 


Ellis et al. (1996) 


-19.20 


-1.09 


2.6 


2.05 


Zucca et al. (1997) 


-19.61 


-1.22 


2.0 


2.2 



a Unit = OMh 3 Mpc~ 3 
b Unit = W 8 hL & Mpc~ 3 
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Table 3. The Baryon Budget 



Component Central Maximum Minimum Grade 3, 



Observed at z ~ 0: 

1 stars in spheroids 0.0026/ifo 1 0.0043/ifo 1 0.0014/ifo 1 A 

2 stars in disks 0.00086/ifo 1 0.00129/ifo 1 0.00051/ifo 1 A- 

3 stars in irregulars 0.000069/ifo 1 0.000116/ifo 1 0.000033/ifo 1 B 

4 neutral atomic gas 0.00033/ifo 1 0.00041/jfo 1 0.00025/ifo 1 A 

5 molecular gas 0.00030/ifo 1 0.00037/ifo 1 0.00023/ifo 1 A- 

6 plasma in clusters 0.0026/ifo 1 - 5 0.0044/if L5 0.0014/if L5 A 
7a warm plasma in groups 0.0056/ifo 1 - 5 0.0115/if L5 0.0029/if L5 B 
7b cool plasma 0.002/ifo 1 0.003/ifo 1 0.0007/ifo 1 C 
7' plasma in groups 0.014/ifo 1 0.030/ifo 1 0.0072/ifo 1 B 



8 sum (at h = 70 and z ~ 0) 0.021 0.041 0.007 



Gas components at z « 3: 

9 damped absorbers 0.0015/ifo 1 0.0027/ifo 1 0.0007/ifo 1 A- 

10 Lyman-a forest clouds 0.04/ifo 1 - 5 0.05/ifo 1 - 5 O.Ol/ifo 1 - 5 B 

11 intercloud gas (Hell) • • • O.Ol/ifo 15 O.OOOl/ifo 1 B 



Abundances of: 

12 deuterium 0.04/if 2 0.054/if 2 0.013/ifo 2 A 

13 helium O.OlO/i™ 2 0.027/1™ 2 A 



14 Nucleosynthesis 0.020/if 2 0.027/if 2 0.013/if 2 



a Confidence of evaluation, from A (robust) to C (highly uncertain) 



